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Abstract 

We study the fluctuation properties and return-time statistics on 
inhomogeneous scale-free networks using packets moving with two dif- 
ferent dynamical rules; random diffusion and locally navigated diffu- 
sive motion with preferred edges. Scaling in the fluctuations occurs 
when the dispersion of a quantity at each node or edge increases like 
the its mean to the power fi. We show that the occurrence of scaling 
in the fluctuations of both the number of packets passing nodes and 
the number flowing along edges is related to preferential behaviour in 
either the topology (in the case of nodes) or in the dynamics (in case 
the of edges). Within our model the absence of any preference leads 
to the absence of scaling, and when scaling occurs it is non-universal; 
for random diffusion the number of packets passing a node scales with 
an exponent [i which increases continuously with increased acquisi- 
tion time window from \i = 1/2 at small windows, to fi = 1 at long 
time windows; In the preferentially navigated diffusive motion, busy 
nodes and edges have exponent fj, = 1, in contrast to less busy parts 
of the network, where an exponent [i = 1/2 is found. Broad distribu- 
tions of the return times at nodes and edges illustrate that the basis 
of the observed scaling is the cooperative behaviour between groups 
of nodes or edges. These conclusions are relevant for a large class of 
diffusive dynamics on networks, including packet transport with local 
navigation rules. 
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1 Introduction 



In transport processes on networks complex dynamical behaviour may be 
caused by the structure of underlying network geometry (for a recent re- 
view see [1] and references therein). The inhomogeneous degree distributions 
found in the scale-free networks may be one cause of this behaviour, how- 
ever, other structural details appear to be essential when dynamics include 
the local navigation of packets [H El E] • Queuing effects at nodes, which is 
essential for high density information packet transport p], El El [U El E] , and 
is usually absent in other transport processes such as charge transport [7] , is 
another cause of collective behaviour in transport dynamics. 

Many dynamically measurable outputs exhibit scaling features, indicat- 
ing a degree of universality that, in turn, can be used to probe a network's 
structural and traffic properties [1J. Distributions of transport times, waiting 
times and noise fluctuations exhibit power-law behaviour in different network 
models. However, it is not an easy task to relate the emergent dynamical 
features to particular structural properties, firstly because of the network's 
power-law inhomogeneity but also because different structural elements, i.e., 
nodes or edges, or higher (hidden) structures, may play a role in the dynam- 
ics. Therefore, a more systematic study of traffic scaling properties, that 
considers structure beyond the node's connectivity, is necessary in order to 
reveal the origin of scaling and relate the observed scaling properties to the 
specific structural and dynamical features on complex networks. 

In view of the long-range correlations in packet streams on inhomogeneous 
networks, analysis of traffic noise, defined by the number of packets processed 
by a node, and traffic flow, the number of packets passing along an edge, 
can give interesting information about traffic conditions and the underlying 
network structure [T]. In recent studies a multi-channel analysis of traffic 
noise on networks [H [9j [TO] reveals that the set of fluctuations {<7j} of the 
traffic time-series {hi(t)} measured at all nodes in the network i = 1,2, ■ ■ ■ N 
obeys the scaling law 

(Ji-Khi^ . (1) 

The exponent /x is often found to be either /i=l/2or/i = l, suggesting a 
super-universal behaviour across different networks [8]. However, recently it 
was found [HI [101 [1] that the exponent \i may depend on traffic conditions 
and the type of measurements taken. Similar results, in which non-universal 
scaling was found to depend on the acquisition time window, where found 
in the analysis of stock market time series [UJ. A unique exponent between 
these two universal values was also found in the analysis of the genome-wide 
time series of the gene expression of yeast [12], the dynamics of which is 
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naturally limited by the cell cycle. 

In this paper we address the question of scaling in the diffusive dynamics 
on networks, by carefully selecting the inhomogeneity of the network and 
the dynamical rule, in order to determine the origin of the scaling and its 
robustness. To demonstrate the importance of these dynamic phenomena 
in real-world networks, we employ a non-trivial model for the transport of 
information packets [TSJ [21 El [I]- For the purpose of this work we use a 
simple scale-free network [14, [15] with a clear inhomogeneity in the node 
connectivity but with low clustering and no edge correlations. We consider 
two types of diffusive motion of the packets on the graph: random diffusion 
and local navigation with preferred edges, as described below. 

In comparison to the models considered in [HI E] , our traffic model is more 
realistic in that: 

• Packets are created at a given rate R and travel to specified destinations 
on network; 

• Packets queue dynamically at nodes; 

• Packets are navigated locally according to a specified algorithm. 

Consequently the travel-times are determined self-consistently rather than 
being fixed as an adjustable external parameter of the model and have a 
broad distribution with a power-law tail, which depends on the network 
structure navigation algorithms [21 [16]. The waiting times of packets in 
queues is another property of our model that is determined self-consistently 
by the dynamics. Depending on the type of queuing discipline employed, the 
waiting time distribution can also have a power-law tail on scale-free networks 
[2]. Throughout this work we use low packet densities in order to keep the 
time series of traffic noise and flow stationary and avoid the effect of large 
queuing times (see [T] for a study of dense traffic on structured networks). 

Our main findings suggest that the presence of a preference in either the 
topology, such as in node connectivity, or in the dynamics, such as by edge- 
preferred navigation, leads to the scaling of fluctuations. When this scaling 
occurs, careful analysis shows that it is non-universal, depending both on the 
acquisition-time and on the importance of the nodes or edges in the transport 
process. These findings of collective dynamical behavior on networks are 
further substantiated with a study of the return-times statistics for nodes 
and edges. 

The organization of this paper is as follows. In section 2 we define the net- 
work and the traffic model and summarize the main features of the transport 
process. In section 3 the results for noise and flow fluctuations are given for 
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various parameters and dynamic rules and in section 4 the origins of scaling 
are discussed and compared with transport on much simpler geometries. In 
section 5 the results for the return-times distributions are obtained for the 
network dynamics studied in section 3, and section 6 gives a short summary 
and discussion of our results. 

2 Traffic of information packets with queuing 
and navigation on network 

2.1 Network structure 

We consider a simple scale-free network [TU [15] with a power-law inho- 
mogeneity in node connectivity, low clustering and no correlations between 
edges. Such a network can be easily grown with a preferential attachment 
rule in which a node with two links is added to the network at each time 
step and each link is connected to a node of degree k with a rate propor- 
tional to k + a. In the emergent structure the degree of the i th added node 
ki ~ {i / N)~ l I ( 1+a \ and hence the connectivity distribution has a power-law 
tail, P(k) ~ k~ T with r = 2 + a, as shown analytically [T5]- A detailed 
analysis of the structure of the networks grown in this way shows that the 
clustering is very low and link correlations are entirely absent, in contrast to 
correlated scale- free networks, grown with the algorithms described in [17j. 
For our simulations of packet diffusion we grow a simple uncorrelated scale- 
free network of N = 1002 nodes and E = 2N edges, which has a connectivity 
distribution with a power-law exponent r « 2.5. 

Once the network is generated, we consider it's structure as fixed, and 
start the transport processes on it. During these dynamics, at each time step 
each node can create a new packet with probability R/N. At creation each 
packet is assigned a destination address, another node on the network where 
it should be delivered. Packets move diffusively through the network, either 
performing a random walk or navigated according to the algorithm described 
below. Once at the destination address packets disappear from the network. 
In this model packets move towards their destination simultaneously, forming 
queues at the nodes on the way. We assume a finite queue buffer of length 
H = 1000 at each node. The first in - first out (FIFO) queuing discipline is 
applied. The maximum length of the queue is important for transport close 
to the congestion state and has an impact on the scaling of the fluctuations 
on nodes [H [TO] . 
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2.2 Random diffusion and edge-preferred navigation 

The motion of packets through the network can either be random or nav- 
igated using some rules, which may affect the role of different nodes and 
edges in the transport process. In the case of random diffusion (RD) a node 
that is processing a packet selects one of its neighbouring nodes at random 
to send the packet to it. Hence the nodes on the network are visited at a 
rate proportional of their connectivity, and hence nodes with a high connec- 
tivity are often visited by moving packets. Implementation of a navigation 
algorithm [T£l [1] will create a bias in the use of nodes and edges on a given 
topology. Here we consider one such navigation algorithm, referred to as the 
CD algorithm, that introduces a preference for less used edges in the traffic 
history [18J. In particular, a node i that is processing a packet at time step 
t selects one of its neighbour nodes, j, that has the minimum value of the 
quantity Sij(t) defined as the product 

Sij (t) = kjF^t) , (2) 

where kj is the degree of the node j and Fij(t) is the number of packets for- 
warded from node i — > j up to the time t. This means that the CD algorithm 
is completely deterministic, involves a search of the nearest neighbourhood 
of the processing node, prefers to send packet to nodes away from the hubs 
and hence introduces a dynamical inhomogeneity to the network transport. 
In this way, from the point of view of nodes, the topological inhomogeneity 
of the network appears to be dynamically reduced. At the same time, an 
unequal use of edges appears, that is measured by the packet flow on them. 
Note that this property of the navigation rules may reduce jamming problems 
at high traffic density, which often occur on scale-free networks in algorithms 
based on the shortest paths [H] • We consider low packet density by keeping a 
low posting rate R, which is much below the jamming rate for this network's 
structure and navigation rules, in order to minimize the potential effects of 
long queues and to retain the stationarity of the traffic time series. 

2.3 Traffic properties 

For the fixed network topology described above and a fixed posting rate of 
packets R, we simulate packet transport both with random diffusion rules and 
with the CD navigation algorithm. The transport properties are measured 
by a number of global and local statistical quantities [131 [21 El U], which 
depend on the diffusion rules. Here we summarize several traffic properties 
which are relevant for further discussion. In Fig. [1] we show how the two 
diffusion rules effect the distributions of travel times of packets and flow 
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along the edges on the scale-free network. Both in the random diffusion and 
in the navigated transport the distribution of travel times exhibits a power- 
law tail, however, both slopes and cut-offs are different. The distribution of 
dynamic flow along edges of the network is shown in Fig. [Tb. For RD the 
dynamical flow is similar to the topological flow (the centrality measure of 
edges) and is dominated by an average value with a width which depends 
on the packet density. In contrast, the CD navigation algorithm induces a 
non-symmetric flow distribution created by the dynamically preferred edges. 

Further differences in traffic with the two diffusion rules are shown in 
the time series of the number of packets processed by individual nodes and 
edges, shown in Fig. [2j A detailed analysis of such time series is given in the 
following section. 
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Figure I: (a) Distribution of travel times of packets and (b) distribution 
of flow along the edges for both random diffusion (RD) and edge-preferred 
navigation (CD). 



3 Scaling of Noise and Flow 

Noise and flow are two local characteristics of the traffic that are determined 
by the number of packets transported in a given time window tk at a node % 
hi(tk), and along a link i — > j and j — > i, fij{%). The index k = 1, 2 • • • K max 
enumerates time windows of length T WIN time steps. 
The flow fij(tk) can be written as 

fij(tk) — Fij(tk) + Fji(tk) ~ [Fij{tk-i) + Fji(tk-i)) (3) 
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Figure 2: Example of time series recorded at (a) a preferred node and (b) at 
an edge with random diffusion rules and time- window Twin = 1000 steps. 

where tp. = kT WIN . During the transport process we record a set of fluctuat- 
ing time-series, one for each of N nodes i = 1, 2, • ■ • N, and similarly 
another set of time-series collected for each of the 2N edges {fij(tk}} on the 
network. In the simulations we fix the creation rates of packets at R = 0.01 
for the RD and R = 0.1 for CD-navigated diffusion, which are well below the 
respective jamming rates on this network structure. 



3.1 Noise fluctuations at nodes with scale-free connectivity 

We determine the dispersion cij and an average < hi > of long time series 
recorded at all nodes {hi(tk)}, i = 1,2, - ■■ N and k = 1, 2, • • • K max . Plots 
of a - , against < hi > are given in Fig. [HK where each point represents one 
node of our scale-free network. The two curves correspond to the results for 
RD dynamics and for the transport with the CD navigation algorithm and 
a fixed acquisition time window Twin = 4000 time steps. As the Fig. [3h 
shows, these plots obey the scaling equation (JTJ with a well defined exponent 
/i, which is different for the two algorithms. For RD the value of fi is clearly 
between the two limits 1/2 and 1 mentioned above, whereas in the case of 
the CD navigation algorithm fi appears to be close to 1. 

By changing the width of the time window in which the data is collected 
we find that the scaling law still holds, but with a different exponent fi. This 
applies for both RD and CD, however, the functional dependencies u(Twin) 
are different. The results are shown in Fig. [Hb- As a rule, if data are acquired 
in a longer time window, the exponent fx tends to be larger. For RD we find 
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Figure 3: (a) Dispersion <7j against average < hi > of time series recorded 
at nodes of the network within a fixed time window Ty/j^ for random (RD) 
and navigated (CD) diffusion, (b) Dependence of the scaling exponent fi on 
the width of the time window. Slope of the line is (3 = 0.11. 

the change from weak dependence at short time windows to steeper increase 
when large windows are applied. For CD the exponent is constantly large for 
large time windows, and changes rapidly when the windows are smaller than 
the typical T WIN « 500 time steps in these simulations. In fact, for CD we 
find that the dispersion between the groups of nodes exhibiting the scaling 
with large fi exponent is gradually reduced with decreasing time window 
and eventually the group becomes so condensed that an exponent cannot be 
defined. Simultaneously another group emerges with the scaling exponent 
close to fi = 1/2. 

3.2 Scaling of flow at preferred edges 

We apply a similar multi-channel analysis to the data on traffic flow fluctua- 
tions. The flow along an edge between nodes i and j is the sum of flow from 
i — > j and from j — > i. In Fig. H]we plot the standard deviations of the flow 
time series against the average flow along that edge < >. Each point 
on the plot represents one edge on the network. As the Fig. shows, in the 
case of CD navigated diffusion, flow fluctuations obtained at a large number 
of edges follow the same scaling pattern as that described in Eq. ([I]) for noise 
fluctuations at nodes. However, in the case of edges, the situation becomes 
different when random diffusion is used: the edges form a dense group on the 
plot, representing almost equal fluctuation properties. This implies that the 
dynamical preference in edges, which is built into the CD navigation rules, 
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introduces both an uneven flow along different links and a different fluctu- 
ation pattern of the flow. The exponent that is measured from the data in 
Fig. UK is close to \i = 1. In the same plot we find another group of edges for 
which the scaling exponent cannot be defined. When the width of the time 
window is reduced the fraction of edges that belongs to the scaling regime 
with \i « 1 is reduced, and a number of edges appear in a new group, for 
which the scaling exponent is close to \l ~ 0.5. The transition between these 
two groups is quite sharp for the deterministic CD navigation (see Fig. HJd). 
We investigate this question further by considering a simplified, probabilistic 
version of the edge-preferred navigation, called D algorithm [18], in which a 
packet moves from node i to its neighbour j with probability 

*4 = 1 " ( 4 ) 



N 



Cijkj 



where CV,- is the adjacency matrix and as before kj is the degree of node 
j. Hence packets are more likely to move towards neighbours with a low 
connectivity. With this navigation rule, we find the flow fluctuations, shown 
in Fig. [5j exhibiting two distinct groups of edges with a smooth transition. 




Figure 4: (a) Dispersion against average < > of time series recorded 
at edges of the network within a fixed time window T WIN = 4000 steps 
for random diffusion RD and CD navigation rules, (b) Same as (a) but for 
different time windows and CD navigation rule. 

Study of the flow fluctuations suggests that the dynamical preference of 
edges, as in our navigation rules, introduces the necessary distinction be- 
tween edges that leads to the scaling behaviour. In contrast to the noise 
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Figure 5: The flow fluctuations against average flow < > for navigated 
diffusion with probabilistic preference for edges pointing to low- connectivity 
nodes. 



fluctuations studied above, in flow fluctuations we did not observe a unique 
scaling exponent with a smooth dependence of the acquisition time window. 
According to their flow fluctuations, all edges appear to be in one of the 
two groups of well defined scaling exponents, either 1/2 or 1. The relative 
population of the groups changes, where the size of the group with the larger 
exponent is increasing with the length of the time window. 

4 Origin of Noise Scaling on Complex Net- 
works 

In the previous section a comparative analysis of the fluctuations of noise and 
flow suggested that in the diffusion processes a preference between network 
elements, that can be achieved in an inhomogeneous network, is necessary for 
the scaling of the type in Eq. CD to occur. In order to further clarify this point, 
we have analysed different geometries including regular and homogeneous 
substrates. 

In Fig. and b we show the fluctuations of packets in random diffusion 
recorded on nodes and edges of a square lattice and of a random graph. 
Clearly, fluctuations of flow at edges on these structures appear to belong 
to one group within a statistical dispersion. However, the record of noise at 
different nodes already shows grouping according to their connectivity: on 
the square lattice one can distinguish between fluctuations in the number 
of packets processed by the four corner nodes, by nodes at the boundaries 
and by the rest of the nodes, whose the connectivity is four. A larger range 
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of groups of nodes is found on a random graph, with connectivities ranging 
from 2 to 8 edges per node, where the larger connectivity nodes are shifted 
to the right part of the plot. In the same spirit, on the scale- free network 
the span of different groups of nodes, with different connectivities is even 
larger, also visible in Fig. [3k, and it is related to the connectivity profile. 
Therefore, the direct relationship between the node connectivity and the 
number of processed packets by that node (its dynamic centrality) occurs 
in the random diffusion. This represents the basis of the observed scaling 
of noise fluctuations on the scale-free networks. In Fig. Eb we show the 
emergent linear dependence between the average noise at a node < hi > and 
its connectivity hi for our scale-free network, which applies for the case of the 
random diffusion algorithm (see also [H]). However, when the diffusion rules 
are changed, as in our navigated diffusion, for example, these relationships 
are altered. We find that in the case of CD navigation the average number 
of packets at nodes of large connectivity plateaus according to the functional 
form < hi >= a — bc ki , with c = 0.6. 

Hence, nodes with different connectivity play a different role in the dif- 
fusion processes on networks. This idea of dissimilarity can be extended to 
the edges of the network when diffusion rules, such as our CD navigation, 
are implemented. Again, we observe distinct groups of edges which have dif- 
ferent flow fluctuations, with more 'important' groups having proportionally 
larger flow and large fluctuations. 




Figure 6: For random diffusion, the dispersion against average of time series 
recorded at nodes and links on (a) regular square lattice and (b) random 
graph. In both figures, the most left group of points represents the data 
collected at (edges links), while the remaining groups are for nodes, (c) 
The average number of packets processed by a node against node degree for 
random (RD) and navigated (CD) diffusion. 

It is an intrinsic property of the diffusive dynamics that different groups 
of nodes or edges develop a cooperative behavior and fall on a line with a well 
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defined profile, as given in Eq. (fT)). Our results for the diffusive dynamics of 
information packets suggest that stationarity of the time series is a necessary 
condition for such cooperativity to occur, whereas the dynamic continuity 
need not be strictly observed (in our model sources and sinks of packets occur, 
which are balanced in the average). Our results show that in topologically 
or dynamically inhomogeneous systems nodes or edges may develop different 
levels of cooperative behavior that leads to non-universal scaling. The non- 
universality is represented by the fact that the scaling exponent depends 
on the dynamics, with either a continuous change between the two limiting 
values 0.5 < p, < 1, or with two exponents defined for different groups of 
edges. A possible origin of these two limiting values of the scaling exponent 
has been discussed in [SI El ED] • 

5 Return-Times to Nodes and Edges 

Another type of dynamic measure collected at individual nodes and edges, 
that depends on the dynamic behaviour of the whole network, is the statistics 
of return-times, or time intervals between the successive events at a given 
node or an edge. In collective dynamical systems such as earthquakes (2TJJ, EH 
[22] . critical sandpiles [231 El], and stock market dynamics [251 EE], a broad 
distribution of return times (sometimes called waiting times or recurrent 
times) is always found, with power-law tails suggesting the occurrence of 
long-range dynamic correlations between the events. In this work we address 
the question of return times to nodes and to edges in order to investigate 
further the nature of collective dynamic behaviour in our model of diffusion 
of packets on a scale-free network. 

Return- Time to Nodes. In the case of random diffusion (random walks) on 
networks the return time distribution has been studied in other parts of the 
theoretical physics literature. In particular, the first return time to the origin 
of a random walker on sparse random graphs, with nodes representing states 
of a system, was considered by Bray and Rodgers [27] as a model of non- 
exponential relaxation in spin glasses and other non-ergodic systems. With 
the help of some heuristic arguments, they arrived at the conclusion that on 
a random graph the long-time behaviour in the diffusion in the phase space is 
dominated by the parts of the network with linear chains (no loops), leading 
to the expression -P(At) ~ exp(— A(p)(Aty^ 3 ), where p is the homogeneous 
connectivity of the random graph and A(p) is known. 

These arguments can be generalized to introduce a power-law distribution 
of connectivities. If the distribution of p behaves like ~ p~ T , and using the 
result in [27] that as p — > 0, P(At) ~ pexp(— 2At/p), then integrating over 
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p leads to P(At) ~ (At)" TA with 



r A = r - 2 . (5) 

Thus, the inhomogeneous connectivity creates a power-law distribution in 
the return times distribution for RD. Recently a more rigorous treatment of 
random walks on scale- free networks was carried out by Noh and Rieger [19] , 
that yielded identical results. The results of our simulations for different 
diffusion processes are shown in Fig. [7h. The return-time distributions in 
different cases studied here seem to have a power-law behaviour before a 
cut-off. (The cut-off can be related to the network size in the case of single 
random walker. Note also a characteristic splitting at small AT with an 
inherent preference for even return times, caused by the lack of clustering 
and the low density of walkers.) In the case of non-interacting random walks, 
i.e., random diffusion without queuing, the results agree, within error bars, 
with the above theoretical prediction. We have the exponent r A = 0.56 ± 
0.03, whereas the distribution of the network's connectivity has a power-law 
exponent r ps 2.5 (see Sec. 2). 

Increasing the traffic density reduces the value of the cut-off, but the 
slope remains practically unchanged. However, when the navigated diffusion 
is considered, both the slope and the cut-off of the distribution are changed. 
In Fig. [7b we show the results for navigated diffusion with the CD algorithm, 
random diffusion with one-depth layer search and RD at low packet density. 
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Figure 7: The all-return time distribution to nodes in the scale-free net- 
work for non-interacting and interacting walks with (a) the random diffusion 
algorithm and (b) for interacting walks with different navigation algorithms. 

Return-Time to Edges. In contrast with the return times to nodes, the 
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situation is entirely different from the point of view of edges on the same 
network. The results are shown in Fig. [HJ We find a pronounced difference 
between the random and navigated diffusion in the tails of the distribu- 
tions. In both cases, however, a unique functional form can be found. For 
larger AT the distributions of the return times to edges can be fitted with 
a q— exponential form, which is often related to non-ergodic behaviour in 
dynamical systems [28J : 



P(At) = B[l-(l-q) 



At 

aT 



V(l-9) 



(6) 



In the case of random diffusion, shown in Fig. [8h, the distribution is very 
close to the exponential form, which corresponds to the q — > 1 limit of Eq. 

In fact, we find q = 1.08 in the case of random diffusion, whereas in the 
case of edge-preferred CD navigation q = 1.33, compatible with non-ergodic 
dynamic behaviour. 
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Figure 8: The all-return time distribution to edges for (a) the random diffu- 
sion and (b) CD-navigated diffusion algorithms. The fit lines are explained 
in the text. 



6 Conclusions and Discussion 

We have performed an analysis of the diffusive dynamics on a scale free 
network with a power-law connectivity distribution but without any other 
form of structural inhomogeneity. Using the transport of information pack- 
ets, where the diffusion rules can be modified in various ways by adjusting 
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the navigation of packets at nodes, we were able to show that several dy- 
namical effects appear to be related to the the microscopic diffusion rules. 
Our approach indicates that these findings will be relevant to more realis- 
tic transport problems on networks, which are very different from simple 
random- walks. 

In particular, we pinpointed the importance of not just topological but 
also dynamic preference to the occurrence of dynamic scaling. We im- 
plemented navigation rules that involve preferences between links, which 
is possible on topologically inhomogeneous scale-free networks. The edge- 
preferential navigation rules appear to dynamically homogenize the network 
(at large connectivity nodes) and yield new dynamical phenomena. We fo- 
cused on two types of scaling behaviour that can be obtained from the point 
of view of individual structural elements, nodes and edges of the network, 
and are potentially related to each other: Scaling of noise and flow fluctu- 
ations, on one hand, and scaling of return-time distributions, on the other. 
While the noise fluctuations and/or return times at nodes have been studied 
extensively on different types of networks [El HE QUI H], our work presents the 
first systematic study of the fluctuations of flow on edges and return times 
to edges, and a comparison with the quantities obtained at nodes within the 
same dynamics. Owing to the different roles that nodes and edges play in 
these diffusion processes on an inhomogeneous network, these comparative 
studies lead to the conclusion that certain types of preferential behaviour in 
either nodes or edges is necessary for the occurrence of scaling. The scaling 
is characterized by the exponents \x and ta and the parameter q. 

In particular, no scaling behavior was found in the fluctuations of flow and 
in the return times to edges for random diffusion. On the other hand, when 
edge-preferred navigation is turned on, both a non-trivial scaling of the flow 
fluctuations and a power-law tail in the return-times to edges is observed. 
The differences in node connectivities, studied in regular, random graph, and 
the scale-free structures, leads to differences in the noise fluctuations both 
in random diffusion and in CD navigated diffusion. Accordingly, we find 
non-trivial distributions of the return times to nodes, that are power-law (up 
to a finite-size cut-off), in agreement with theoretical predictions [271 [T§] . 
The exponents, however, depend on the navigation rules. Only in the case of 
random diffusion do the return times scale with a power which is determined 
by the network's connectivity distribution. 

Furthermore, we show that the scaling properties of noise fluctuations 
are non-universal, with both the topological inhomogeneity but also details 
of the dynamics playing a role in the emergent scaling behaviour. In addition 
to a well pronounced dependence on the width of the time window, that was 
also observed in other studies [9], (HE EH [H] , we found that navigation rules 
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effect the universality of noise fluctuations. Again, we found a qualitative 
difference between the scaling properties at nodes and at edges. A unique 
scaling exponent h{Twin) is always found in case of noise fluctuations at 
nodes, whereas in our navigated diffusion the flow fluctuations at edges are 
bi-universal. 

Finally, we have demonstrated that by looking beyond the random walk 
dynamics on sparse topologies, one may find a number of new dynamical 
phenomena, that are both interesting from the point of theory, but also 
practically important for many real transport processes on networks. 
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